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Abstract 



Representation of a D-dimensional fermion determinant as a path integrai 
^ I of exponent of a (D + l)-dimensional Hermitean bosonic action is constructed. 

o 

<N 1 Introduction 

in . _ 

Q^ . In the recent paper we proposed a method which allows to present a four di- 

^ I mensional fermion determinant as a path integrai of exponent of a five dimensionai 

constrained bosonic action. This method provides a path integrai bosonization of 
fermion models in dimensions D >2. The problem of bosonization of fermion mod- 
^ ■ els in dimensions D > 2 was studied previously by several authors (see for example 

IH' @) 01) 05 0)' however no quite satisfactory solution was known. Bosonic rep- 
^ i resentation is particulary important for lattice gauge theories like lattice QCD as it 

I allows to avoid complications related to integration over grassmanian variables. One 

way to handle the problem was proposed in papers 0, [H, where the algorithm for 
approximate inversion of the QCD fermion determinant was formulated replacing 
it by an ifìnite series of bosonic determinants. In our paper the exact effective 
bosonic action was explicitely constructed and it was proven that when the model is 
considered on the finite Euclidean lattice the path integrai is convergent. Any gauge 
invariant observable in QCD may be calculated in this approach. However numerica! 
simulations by Monte-Carlo method are not straightforward as the effective bosonic 
action is not Hermitean. 

In the present paper I propose an alternative procedure which provides the rep- 
resentation of a D-dimensional fermion determinant in terms of a path integrai of 
[D + l)-dimensional Hermitean bosonic action. This method seems to be more ap- 
propriate for Monte-Carlo simulations. The construction can be inverted providing 
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the representation of bosonic determinant in terms of a path integrai of fermionic 
effective action. 

In the next section I illustrate the idea by the model on an infinite lattice. In 
the last section a corresponding construction for a finite lattice is presented. 



2 Infinite lattice model 

To illustrate the idea I firstly present a formai construction for the infinite lattice 
theory. In this formai proof I ignore the problem of convergence of the path integrals. 
In the next section a rigorous proof will be given for the case of a finite lattice. 
Let det(i?) be the determinant of some positive bounded Hermitean operator B. 

det{B) = I exp{-a^^[^„(x)S'">„(x)]}#d^ (1) 

Here il)m{x) is a fermion field, and m is a collective index which numerates spinorial, 
colour components e.t.c. These fields are defined on a D-dimensional Euclidean 
lattice with the lattice spacing a. Our goal is to construct a representation for the 
determinant of the operator S as a path integrai of the exponent of a Hermitean 
bosonic action. 

Let US introduce bosonic fields (t){x, t) which depend on one extra space coordinate 
t and have the same spinorial and internai structure as the fields ■0. The extra 
coordinate t is defined on the infinite one dimensionai lattice with the lattice spacing 
b: t = nb, —oo < n < +oo. The fields are described by a locai Hermitean action 

S = a''b £ ^[. >:+i(^)e'''Vn(x) - 0:(x)e-^V„+i(x) 

n=-oo X 2Ò 

We shall prove the following identity 

det(S) = lim / explS + i{(t)l{x)x{,x) + x*(a;)0n(a;))]}#*(i0„(ix*cix (3) 

Here x(a;) are Z^-dimensional bosonic fields which play the role of Lagrange multi- 
pliers imposing the constraints 

^0„(x) = ^C,(^) = O. (4) 

n n 

R.h.s. of eq. (0) can be written in terms of eigenvectors of the operator B, B^ being 
the corresponding eigenvalues 

r +°° rh*°' piB^bia _ J^*a -iB°'bia 

/ = Ihn/ exp{6 E E r-"^ ^" + 



-oo a 



To calculate the integrai (|^) we make the following change of variables: 

C^exp{z5"nò}C, ^ exp{-2S-nò}C* (6) 
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After this transformation the integrai (§) written in terms of Fourier components 

+00 

<Pk = b (pnexp{-iknb} (7) 

n=— 00 

looks as follows 

/ = lim / exp{-{2n)-^J2 I ! dk[\^k? siYi{kb)b-^ + 

+^{4>Tx'' + x^'^rMk - B^)]}d4>Td4>tdx*''dx'' (8) 
Integrating over (/)„ one gets 

I = limZ J exp{-^Y.Xa[MBab)b-']-'Xa}dx*JXa (9) 
where Z is the Ba independent (infinite) Constant 

Z = f exp{-^Y. r dk\Mk)\'sm{kb)b-'}d4>l{k)dMk) (10) 

This Constant may be ehminated by a proper normahzation. 

If the eigenvalues are hmited we can choose b sufficiently small, so that 
Bab << 1. Therefore we can replace sm{Bab) by Bab, and integrating over Xa to 
get 

I = llB^ = det{B) (11) 

a 

The eq.(D is proven. The determinant of the positive bounded operator B is pre- 
sented as the path integrai of the purely bosonic Hermitean action. 

However the derivation in this section was formai as we did not study the problem 
of convergence. In the next section I shall present an analogous construction for a 
finite lattice, where we are able to prove the convergence of ali the integrals and 
make a representation like eq.(|^) quite rigorous. 



3 Finite lattice models. Lattice QCD 

Having in mind applications to QCD in this section we consider two fermion fiavours 
interacting vectorially with the Yang- Mills fields 1/^. The reason to consider two de- 
generate fiavours is the positivity of the square of the gauge covariant Dirac operator. 
It allows to prove the convergence of the path integrai for the model defined on the 
finite lattice. Formally the construction goes for any number of fiavours, but a rig- 
orous proof will be given for the case of even number of fiavours. Firstly we reduce 
the problem to calculation of the determinant of Hermitean operator by using the 
identity 

det(I) + m) = det[75(Ì) + m)], (12) 
D = ^^.{D; - D,) (13) 
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In eq.(|T3D is the lattice covariant derivative 

D^^{x) = ^[U^{x)^{x + a^) - iPix)] (14) 
is a lattice gauge field. 

It is convinient to present the square of fermion determinant in the following 
forni 

/ expla"^ J2 H Ì'n{x)iD + m)ìpnix)}di>dìp = det[75(D + m)]^ = (15) 

n=l X 

expjg^ ìj){x){D'^ — m^)ìj)[x)}dil)d'ìjj 

X 

The operator —D^ + is the square of Hermitean operator therefore ali it's eigen- 
values are positive . 

We again introduce five dimensionai bosonic fìelds 0(x, t)with the same spinorial 
and internai structure as iJj{x). The spatial components x are defined as above. The 
fifth component t to be defined on the one dimensionai lattice of the length L with 
the lattice spacing h: 

L = 2Nb, -N <n<N (16) 
We choose h « a. The free boundary conditions in t are assumed: 

0n = 0, n< -N, n> N (17) 

The folowing identity is valid 

I exp{a^^^/^(x)(l)2-m2)^/>(x)}#dV'= (18) 

X 

r ^ 

= \im exp{a% J2 E[(2&')"'(C+i(a;) exp{z75(Z^ + m)ò}0„(x) + 

L^oo,b-.OJ n=-N+l X 

+0* (x) exp{-275(£) + m)6}0„+i(x) - 20*0^) + 
+-^{(p*n{x)xix) + x*ix)(/)n{x))]}d<j)ld<j)ndx*dx- 

R.h.s. of eq. (|1^) can be written in terms of eigenvectors of the operator ^^{D + m) , 
which are denoted as Da- 

To check the convergence of the integrai ( |r^) we shall write it in terms of Fourier 
components. For this purpose it is convinient to replace the model on the lattice 
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with 2N points and free boundary conditions by the equivalent model on the lattice 
with 2N + 1 points, periodic boundary conditions and the additional constraint 



-N — (P-N 



(20) 



These constraints may be imposed by adding to the action the corresponding terms 
with Lagrange multipliers A. Then the integrai ([19D may be written in terms of 
Fourier components 

7 2mkn 

0fe = 2^</>nexp{— — — } (21) 



.N 2ÌV + 1- 



as foUows 

^ 2Tlk 

I0fcr[cos(- 

a k=-N 



^ = JZ^oJ l<^^r[cos(^^-D"6)-l](26^)-i+ (22) 



1 ~ ~ ~ 2Ì7rkN 

+z6(L)-2 (^-X" + ;^*>-) + tb{2N + 1)-' ^(0-exp{^^^}A„+ 

—2inkN 

+A>gexp{ ^^^^ })]}d(f)Uk)dMk)dx*adXadXldXa 

The representation (p2|) makes the convergence of the integrai obvious. 

Now we come back to coordinate space and to calculate the integrai (0) we 
make the foUowing change of variables: 

< ^ exp{2D"nò}(/)^, (j)^* ^ exp{-iD"nb}(j)^* (23) 

A° ^ exp{iD"A^6}A", A*" ^ exp{-iD"iVÒ}A*° 

Then one gets 

/= hm /exp{(,E|«l«±^£|±i^^+ (24) 

L^oo.fe^Oj '^AT 2b 

n=—N 

+z(L)-^(C"e^''"'^V + X*"e-^^"'^''C)] + 
+^(&A" + A*°Cjv)}^^0r#n^^X*"rfx°t^A*"dA° 

Now the quadratic form in the exponent does not depend on Da and therefore the 
corresponding determinant is a trivial Constant. So we can calculate the integrai 
by finding the stationary point of the exponent, which is defined by the foUowing 
equations: 

b-\4>Zi + C-1 - 20r) + 2L-^x*"e-^^""^ = 0, n^-N 

r2(C+i + C-i-2C) + ^^"^X°e^''""'' = 0, n^-N (25) 
h-\<ì)*^j^+^ - 20!"^) + iL-'2x*''é^-^^ + = 
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For small b these equations can be approximated by the differential equations: 

4>*'' + iL-^X*"e-'^'^' = (26) 
0" + iL-'^x"e'^"' = 

0(§) = 0(-§) = o, 0*^) = = 



The solution of these eq.s is 



^/L{D 



a\2 



X*" .2t . .D-L D-L 
:[— sm(— — ) + i cos(— — )] 



a\2 



X" ,2t . . D-L 

sm( — - — )— zcos( — - — )J 



Substituting this solution to the eq.( |4D we get 

/ = ^im j exp{-^^^ + O(L-i)} = \{{D-f = det{-D' + m') (28) 

The eq.(|T8|) is proven. It is clear from eqs.(^,^), that in the limit L ^ oo the 
dependence on the boundary conditions disappear. 

For practical calculations it may be useful to linearize the eq.(|T8|). It can be done 
in a straightforward way using the fact that the operator 75(-D + m) is bounded 0: 

\\j5{D + m)\\ < 8a-^ + m (29) 

. Therefore for any a the product -D"ò << 1. The exponent exp{ibDa} may be 
replaced by 1 + ibDa — \b'^{Da)'^ and linear ized version of the eq.([T8|) looks as follows 

I exp{a^^ V'(a;)(^^ - m^)^l){x)}di)dil) = (30) 

X 

N 

= lun exp{a% ^ 5:[(2ò2)-i(C+i(x)0„(x) + 0:(x)0„+i(x) - 20:0„) + 

n=-N+l X 

+i[(Pl_^^{x)-f5{D + m)0„(x) - C(x)75(D + m)<Pn+i{x)][2b]-' + 

1 i 
+ 2^niD^ - m^)(t)n + -^{(t)*n{x)x{,x) + X* (a^) 0„ (x) )] }c/(/);rf(/)„dx*o?X- 

One can get rid off the last terni in the exponent by solving explicitely the constraints 

N 

E 0n(x)=O (31) 

n=-N+l 

Substituting the solution of this equation to the integrai (|^) we get a representation 
for the determinant of the square of the gauge covariant Dirac operator as the path 
integrai of the exponent of the purely bosonic Hermitean action. 
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4 Discussion 



Comparing the results presented above with the previous ones ^ we note that 
the new bosonic representation for gauge covariant Dirac operator is based on the 
Hermitean action and therefore is better suited for numerical simulations. We also 
got rid off the extra parameter A which was present in our first paper. One can hope 
that the construction given in this paper will allow to use directly the Monte-Carlo 
method for calculation of fermion determinant. 
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